Accepted to 2011 March 10 ApJ 

Preprint typeset using I^FfjjK style cmulateapj v. 2/16/10 



THE BINARY FRACTION OF LOW MASS WHITE DWARFS 
Justin M. Brown 

Franklin and Marshall College 
415 Harrisburg Avenue, Lancaster, PA 17604 USA 



o 

< 

o 

Oh- 
i < 

O ■ 



> 

OS 

in 



13 



MUKREMIN KlLIC, WARREN R. BROWN, AND SCOTT J. KENYON 
Smithsonian Astrophysical Observatory 
60 Garden Street, Cambridge, MA 02138 USA 
Draft version January 28, 2011 

ABSTRACT 

We describe spectroscopic observations of 21 low- mass (<0.45 M Q ) white dwarfs (WDs) from the 
Palomar-Green Survey obtained over four years. We use both radial velocities and infrared photometry 
to identify binary systems, and find that the fraction of single, low-mass WDs is <30%. We discuss 
the potential formation channels for these single stars including binary mergers of lower-mass objects. 
However, binary mergers are not likely to explain the observed number of single low-mass WDs. Thus 
additional formation channels, such as enhanced mass loss due to winds or interactions with substellar 
companions, are likely. 

Subject headings: White Dwarfs — stars: Low-Mass — stars: Evolution 



1. INTRODUCTION 

Around 10% of white dwarfs (WDs) in the solar neigh- 
borho od are low-mass (M < 0.45 M©) helium-core ob- 
jects (jLiebert et al.l 120051) . The oldest globular clus- 
ters in the Galactic halo are cu rrently producing «0.53 
Mq WDs (|Kalirai et al.l 12009( 1. Hence, the presence 
of low-mass helium-core WDs in the Galactic disk re- 
quires a formation mechanism other than the normal 
single star evolution. The most commonly accepted for- 
mation channel for low-mass WDs is through enhanced 
mass loss from post-main-sequence stars in interacting 
binary systems. Such stars can lose their outer en- 
velopes without reaching the asymptotic giant branch 
and without ever igniting helium, thus ending up as 
helium-core WDs. The large number of low-mass WDs 
found in binary systems support this formation scenario 
(|Rebassa-Mansergas et al.l 1201 ID . However not all low- 
mass WDs are found in binary systems ( Marsh et al.l 
119951: iMaxted et alJ|2000t iNapiwotzki et al.l 120041) . 

iNelemans fc Taurisl (|1998f ) propose a scenario for the 
formation of single low-mass WDs in which binary inter- 
actions of post-main-sequence stars with close-in mas- 
sive planets can expel the stellar envelope. The plan- 
ets may or may not survive the common envelope evo- 
lution. In the latter case, there is no way of knowing 
whether low-ma ss WDs are produced thro ugh this for- 
mation channel. INelemans fc Taurisl (|1998l ) predict that 
the minimum companion mass required to expel the en- 
velope of the progenitor 1 M w star i s 21 M j. Based on 
mid- infrared photometry. iKilic et al.l ((2010T ) limit any un- 
detected companions to M < 40Mj around eight single 
low-mass WD systems. Therefore, a well tuned common- 
envelope phase scenario involving 20-40 Mj companions 
is required to explain the radial velocity and infrared 
observations of the apparently single low-mass WDs. 
Even though this scenario seems unlikely, it cannot be 
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ruled out due to the caveats in our understanding of the 
common-envelope phase evolution. 

An alternative scenario fo r the formatio n of single low- 
mass WDs was proposed by Hansen (2005) to explain the 
color-magnitude distribution of the WDs in the metal- 
rich cluster NGC 6791. Metal-rich stars have increased 
opacities and thus increased mass loss winds; extrapolat- 
ing theoretical mass loss models to [Fe/H]=+0.4 suggests 
that mass l oss can approa ch 0.4-0.5 M© on the red gi- 
ant branch (lCatelanll2000( l. Thus a single low-mass WD 
may result from a single metal-rich star bec ause of severe 
mass loss winds on the red giant branch. iKalirai et al.l 
(2007) observe that most of the brightest WDs in NGC 
6791 are indeed low-mass WDs, and they suggest that 
at least 40% of the stars in that cluster have lost enough 
mass on the red giant branch to end up as He-core WDs. 
The relatively large fraction of low-mass WDs in NGC 
6791 compared to the 10% fraction in the field WD sam- 
ple likely points to the dependency of stellar evolution on 
metallicity. Thus a mechanism exists to produce single 
low-mass WDs in metal-rich environments. 

In the mass-loss wind scenario, the formation rate of 
single ~0.4 helium WDs is dominated by the old stel- 
lar population; intermediate-mass metal-rich stars, even 
with enhanced mass loss winds, have enough mass to ig- 
nite helium in their cores and eventu ally form ordinary 
carbon-oxygen WDs. According to the lReid et all (2007) 
analysis of stars within 30 pc, 3% of solar neighborhood 
stars with >8 Gyr ages have [Fe/H]>+0.3. This is com- 
parable to the fraction of single l ow-mass WDs in the 
PG survey WD population, which IKilic et al.l (|2007f) es- 
timate is 4%. Thus single metal-rich field stars provide 
another possible origin for single low-mass WDs. 

All three scenarios mentioned above - mass-loss due 
to binary companions, planets, and winds - likely con- 
tribute to the population of low-mass WDs, the question 
is by how much. The binary formation channel is gener- 
ally considered to be the dominant mechanism, however 
the remaining two channels are also likely to contribute. 
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M, T, 



effj 



and V are from ILiebert et al.l ( f2005T) ; J 



from Cutri et al. ( 200J|) , except for those objects (*) with our own 
photometry. 

Understanding the frequency of binary /single low- mass 
WDs is necessary to constrain the formation models. 

In this paper, we measure the binary frequency of the 
30 lo w-mass WDs found in the Palomar Green (PG) Sur- 
vey ([Green et al.l ll986). We provide 720 radial velocity 
observations for the 21 previously unstudied low-mass 
WDs, and investigate the near-infrared color excess of 
the whole sample. We find that at least 70% of the PG 
WDs with < 0.45 M are binaries, and compare this ob- 
served binary fraction and orbital period distribution to 
other samples of WDs and to low-mass WD formation 
models. 

2. OBSERVATIONS AND TECHNIQUES 

Our sample of thirty < 0.45 Mq low-mass WDs is 
taken from the PG Survey, a comp rehensive survey o f 
blue stellar objects with B < 16.1 (iGreen et al.l I1986T ). 
The PG Survey covers 10,714 s quare degrees with an es- 
timated completeness of 84%. ILiebert et al.l (|2005l ) fit 
stellar atmosphere models to the 348 DA WDs in the 
PG survey and identify 30 WDs with mass less than 0.45 
Mq . This is the sample of WDs studied h ere. The mass 
and e ffective temperature determined by ILiebert et al.1 
:2(I(I5; are listed in Table [U 

2.1. Photometry 

We take ]/-band photometry from ILiebert et ail (2005) 
and near-infrared pho tometry from the Two Micron All 
Sky Survey (2MASS) (|Cutri et al.ll200l . In addition, for 
nine objects we obtain deeper near-infrared photometry 



using the Peters Automated Infrared Imaging Telescope 
(PAIRITEL) - the old 2MASS nort h telescope operate d 
with the original 2MASS camera (|Bloom et al.l 12001 . 
The photometry of the 30 low-mass WDs is tabulated 
in Tableffl 

2.2. Spectroscopy 

We obtain spectroscopy for the 21 low-mass WDs with 
no previously known radial velocity variability. The 
other 9 WDs are known b i nary systems published else- 
wher e dMarsh et all 119951: [Marshl 119951: iHolberg et al.l 
19951: lOrosz et all 119991 : (Morales-Rueda et al.l 120051 : 
Nelemans et al . 2005|). To verify our data reduction and 
analysis procedures, we re-observed one of the known 
binaries, PG 2331+290, and confirmed its 4 hr orbital 
period. 

We obtained 720 spectra between October 2007 

and April 20 10 using the FAST spectrograph 

(jFabricant et al.l 119981 ) on the Fred Lawrence Whipple 
Observatory 1.5m telescope. The spectrograph was op- 
erated with a 1.5" slit and a 600 line mm^ 1 diffraction 
grating. This spectrograph set-up provides a wavelength 
range of 3500 to 5500 A at a spectral resolution of 1.7 
A. All observations were obtained with a comparison 
lamp exposure following the science target exposure. 

We process the spe c tra u sing IRAIt] following the 
guidelines in iMassevi (|1997[ ). We determine radial 
velocities using the c ross-correlation package RVSAO 
(|Kurtz fc Minkl Il998f) with the following procedure. 
First, we measure preliminary velocities by cross- 
correlating with a high signal-to-noise WD template of 
known velocity. Second, we shift each object's spec- 
tra to rest frame and sum them together to create a 
high signal-to- noise template of each object. Finally, we 
cross-correlate the individual spectra with their respec- 
tive template to obtain radial velocities with the highest 
possible precision. The average uncertainty of our radial 
velocities is ±16 km s _1 . 

We verify the accuracy of our velocities by measuring 
two night sky emission lines, Hg 4358.335 A and 5460.750 
A. We find that the night sky lines have a km s _1 
mean and a 10 km s _1 dispersion in the data obtained 
in 2009 and 2010. However, the night sky lines exhibit a 
larger dispersion and non-zero means up to ~10 km s _1 
in the data obtained prior to 2009. We account for this 
systematic error by adding ±10 km s~ x in quadrature 
to the radial velocity errors for spectra obtained prior to 
2009. 

3. RESULTS 

We identify binary systems among our low-mass WD 
sample in two ways: by radial velocity variability, and by 
infrared color excess. We cannot detect all binary com- 
panions, of course. Radial velocity variability is sensitive 
to companions of any type, but only for small orbital 
separations and non-zero inclinations. Infrared photom- 
etry is sensitive to companions of any orbital separation 
and inclination, but only for main sequence (M dwarf) 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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stars that produce a significant infrared color excess. De- 
spite these limitations we find that at least 70% of our 
WDs have binary companions, and for these systems we 
estimate the mass and nature of the companions. 

3.1. Spectroscopic Binaries 

We begin by fitting orbits to our radial velocity data, 
and then testing whether the observed velocity variability 
is statistically significant. We search for orbital periods 
in our radia l velocity data using a Lomb periodogram 
(|Presd fl994l) . Period aliases are present for all of our 
WDs, thus we select the peri od that minimizes y 2 fo r 
a circular orbital fit following iKenvon fc Garcial (1986). 
Figure [T] (see also Figure [3]) plots the radial velocities 
phased to the best-fit periods for the objects with well- 
determined orbits. 

We use an F-test to calculate the significance level at 
which the orbital fits have a smaller variance than a con- 
stant velocity fit. For a significance threshold of 0.01, 
we find 7 WDs with significant orbital fits. The orbital 
elements for the 7 WD s are p resented in Table El As dis- 
cussed in lBrown et al.l ((2010), uncertainties in the orbital 
elements are derived from the covariance matrix and \ 2 ■ 

Curiously, of the WDs for which we cannot fit orbits, 
two objects (PG 0846+249 and PG 1320+645) have sig- 
nificant velocity variability (listed in Table El as having 
"high dispersion"). Although we require more data to 
constrain a period for these two WDs, we consider them 
probable binary systems. 

Table El summarizes the orbital properties for the entire 
sample of 30 low-mass WDs. For our 7 newly discovered 
binaries, we present the orbital period P, velocity semi- 
amplitude K, systemic velocity 7, time of spectroscopic 
conjunction To, the F-test significance, the mass function 
MF, the minimum companion mass M2, and the maxi- 
mum merger time r. For the 14 WDs with no orbital 
fit, we present the systemic velocities and, in place of the 
semi-amplitude, the standard deviation of the velocities. 
Finally, for the 9 previously known binaries, we present 
the published periods and mass functions as noted in 
Table El 

3.1.1. Orbital Properties 

Sixteen of the 30 low-mass WDs in our sample (7 from 
our observations and 9 from the literature) are binaries 
with well-measured orbits. For these systems we calcu- 
late the companion masses, merger times, and likelihood 
to produce Type la supernovae (M tot > 1.4 M Q ). 

The orbital mass function relates the observed period 
P and velocity semi-amplitude K to the mass of the two 
components, M\ and M2, and the inclination i: 



PK 3 (M 2 smif 



2ttG (Afi + M 2 ) 2 ' 

Because the mass of th e WD, M-\ , is known 
atmosphere model fits (jLiebert et al J 120051) . 
the mass function to solve for the mass of 
companion, M2 sini. An edge-on orbit with i 
the minimum possible companion mass and 

in Table El 

Our WD binaries are potential candidates 
A short-period binary loses energy through ; 
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Fig. 1. — Radial velocities for 6 WDs with well-determined orbits, 
phased to the best-fit periods (see Table 0. 



wave radiation, leadin g to an inward spiral of th e ob- 
jects until they merge (jPostnov fc Yungelsonll2005l). We 
calcul ate the merger times following lLandau fc Lifshitzl 

cin), 



T = 



(M!+M 2 )3 
MiM 2 



?5 x 10^ 2 Gyr, 



(2) 



where the masses are in Mq, and the period is in hours. 
We use the minimum companion masses to calculate the 
maximum possible merger times, presented in Table El 

Two of the WD+WD binaries, PG 1101+364 and PG 
2331+290, will merge within a Hubble time. If the com- 
bined mass of the system exceeds the Chandrasckhar 
limit, the merger may result in a Type la Supernova. The 
cumulative probability function of orbi tal inclination for 
a ran dom stellar sample is 1 — cosi (jlmbert fc Prevotl 
1998). Using this relation, we calculate the probability 
of the merging binary systems becoming Type la Super- 
novae. None of these systems have probabilities greater 
than 2% of becoming Type la Supernovae. Two ad- 
ditional systems in our survey, PG 1224+309 and PG 
1458+172 (discussed below) are WD+MS binaries with 
merger times shorter than a Hubble time. These two sys- 
tems will evolve into cataclysmic binaries within several 
Gyr. 

3.2. Photometric Binaries 

We now use infrared photometry to search for binary 
companions to the low mass WDs. We begin by com- 
paring the observed (V — J) color of eac h WD with the 
expec ted (V—J) color from the models of lBergeron et all 
(1995) for the appropriate WD mass and effective tem- 
perature. If there is a (V—J) color excess with more than 
2cr significance, we classify the WD as having a probable 
(M-dwarf) companion. 
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Fig. 2. — The spectral energy distributions for four WDs with infrared excess. SDSS and 2MASS photometry are plotted with red dots. 
Spectroscopy is drawn in black. The best-fit WD (dashed red line) plus M dwarf model is drawn with the solid red line. 



Five of the 30 WDs have significant infrared color 
excess, and are listed in Table [5J A sixth WD (PG 
0237+242) with possible color excess was shown by 
iKilic et ail (|20100 to have n o mid- infrared signature 
of a companion. IKilic et all pOlOh also studied PG 
2257+162, one of our five WDs with color excess, and 
found clear evidence for a 3300 K (« M4 spectral type) 
companion in the mid-infrared Spitzer photometry. 

Figure [2] plots the spectral energy distribution for the 
other four WDs with color excess. All four objects also 
have ugriz photometry available from the Sloan Digi- 
tal Sky Survey. In Figure [5] we compare the optical 
and near-infrared photometry (points) to the observed 
SDSS spectra (black line) and to theoretical WD plus M- 
dwarf models (red lines) . All four WDs are well modeled 
with main seq uence M-dwar f companions with spectral 
types M3-M5 (|PickleslH99l ). Based on known M3-M5 
dwarfs in eclipsing binaries, we estimate that M3, M4, 
and M5 dwa rfs have masses of 28, . 22, and 0.19 M n 
respectively dCakirli fc IbanogliJl2010t llrwin et"alll2010t 
iMorales et afl l2009). The formal precision of these mass 
estimates is ±0.05 M Q . 

3.3. Comparison of Photometry and Spectroscopy 

Knowing the masses and the velocity amplitude (or 
limits to the velocity amplitude) of the binary systems, 
we can place limits on the orbital periods and inclina- 
tions using the mass function (Equation [1]). We separate 
period and inclination by assuming, for purposes of dis- 
cussion, the mean inclination angle of a random stellar 
sample, i = 60°, and the mean orbital period observed in 
our sample, P ~ 1 day. Table [5] summarizes the results. 
The columns in Table [5] are M^phot, the photometric 
mass estimate, P Foi i= g ° > the orbital period given the 
observed velocity amplitude and a 60° inclination, and 



TABLE 3 
Photometric Binaries 



Object 


A/ 2 ,phot (M Q ) 


P For l = 60° ( da y S ) 


'■For P = lday (°) 


PG 1224+309 


0.22 






PG 1252+378 


0.19 


5.3 


30° 


PG 1415+133 


0.22 


20.4 


18° 


PG 1458+172 


0.28 






PG 2257+162 


0.22 


8.0 


26° 



ipoi p=iday, the orbital inclination given the observed ve- 
locity amplitude and a 1 day period. 

Three systems with color excess have no detected radial 
velocity variability. The implied orbital parameters for 
these three systems, given the photometric companion 
mass estimates, are consistent with the radial velocity 
observations. Given the 20-30 km s _1 upper limits to the 
velocity semi-amplitudes, PG 1252+378, PG 1415+133, 
and PG 2257+162 could be binaries with either week- 
long orbital periods (assuming i = 60°) or relatively pole- 
on i < 30° inclinations (assuming P=l day). We would 
not detect any of these systems with our radial velocity 
data. 

Two systems with color excess, PG 1224+309 and 
PG 1458+172, have significant radial velocity variabil- 
ity. Both systems show hydrogen and magnesium emis- 
sion lines from the irradiated face of the secondary star. 
IQrosz et al.l (|1999D used radial velocity observations of 
PG 1224+309 to constrain the mass of the companion to 
0.28 ± 0.05 Mq. This mass estimate is consistent with 
our photometric mass estimate within the errors. Given 
the 6.2 hr period of this binary system, PG 1224+309 
will evolve into a cataclysmic variable system within a 
Hubble time. 
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Fig. 3.— (Top) Typical PG 1458+172 spectrum, unsmoothed, 
with H and Mg lines seen in emission. (Bottom) Observations 
phased to the best-fit 3.968 hr period for both absorption line ve- 
locities (squares) and emission line velocities (x's). 

3.4. PG 1458+172 

PG 1458+172 is a new pre-cataclysmic binary discov- 
ered in our survey. The optical spectrum of this ob- 
ject shows Balmer lines both in absorption and emission 
and a Mg I triplet (A 5167.321 A , 5172.684 A , and 
5183.604 A) in emission as seen in Figure [H Strong Ha 
emission is also evident in the SDSS spectroscopy. The 
radial velocity variation of the emission lines is about 
190° out of phase with the Balmer absorption lines from 
the WD primary. We use the H and Mg emission lines 
to fit an independent orbit to the secondary, and list 
the orbital parameters in Table [5] under PG 1458 (Mg). 
Both the primary and secondary have best-fit periods of 
0.16532 ±0.00033 days. 

The ratio of the velocity a mplitudes, combined with 
the 0.41 Mq WD mass from (|Liebert et all 120051 ). sug- 
gests that the secondary has a mass of 0.55 M Q . How- 



ever, the center of light of the emission features is dif- 
ferent than the center of mass of the companion. If the 
emission features come from the heated face of the com- 
panion, for example, we would expect a shift in the veloc- 
ity se mi-amplitude of the emission features (|Orosz et al.1 
1999). Fortunately, we detect the light of the compan- 
ion in the infrared (see Figure [5]) which provides a more 
reliable companion mass estimate. 

Figure [2] shows that the infrared excess around 
PG 1458+172 is consistent with an M3 (0.28 M©) dwarf 
companion. This is significantly lower than predicted by 
the mass function and the spect r oscop ic mass measure- 
ment for the WD. iLiebert et al.l (|2005l ). unaware of the 
Balmer emission features, likely over-estimated the WD 
mass, however. Higher signal-to-noise ratio spectroscopy 
covering the Ha emission line is needed to properly char- 
acterize the orbits and masses of the individual compo- 
nents in this pre-cataclysmic binary. 

3.5. Binary Fraction and Completeness 

Of the 21 low- mass WDs we observed from the 
Palomar-Green Survey, seven of these have clear periodic 
radial velocities, two have large enough dispersion to be 
binary systems (without our being able to constrain a pe- 
riod) , and three of those without strong velocity variation 
have infrared excess, indicating a companion. Including 
the nine objects from the literature that comprise the 
30 low-mass WD sample of the Palomar-Green Survey, 
21 show significant evidence for a companion. Thus the 
binary fraction of our ~0.4 Mq WDs is at least 70%. 

Some binaries must escape our detection, however. We 
are not sensitive to binary systems with intrinsically faint 
companions (such as brown dwarfs and WDs) that have 
pole-on orientations, or large orbital separations. For a 
sample of 30 binary systems with random orbital incli- 
nations, there should be 4 pole-on systems with i < 30° 
for which we would not detect radial velocity variability. 
Given our discovery rates, one of these hypothetical pole- 
on systems should be identified by infrared photometry. 
Thus there may be an additional 3 (10%) more binaries 
in our sample, undetected because they are pole-on sys- 
tems. Put another way, 20% - 30% of our low-mass WDs 
may exist in single systems. 

Additional observations su ggest that some of th e low- 
mass WDs are indeed single. iMaxted et al.l (|2000l) have 
7-15 radial velocity measurements for five of the possi- 
ble single systems (PG 0132+254, 0808+595, 1229-013, 
1614+136, and 2226+061) and do not detect signif- 
icant radial velocity variations for any of these ob- 
jects. Equally im portant, mid-infrared photometry by 
iKilic et all (|2010D rules out the presence of both stellar 
and massive brown dwarf companions in six of the ap- 
parently single systems in our sample. We conclude that 
possibly six (20%), but no more than nine (30%), of our 
30 targets are single low-mass WDs. 

4. DISCUSSION 

4.1. Comparison to Population Synthesis Models 

The observed fraction of WD binary systems is related 
to WD mass. WDs with 0.5-0.7 M Q a nd 0.7-1.1 M ff 
observ ed in the volume- limited sample of iHolberg et al.l 
(2008) have binary fractions of at least 32% and 6%, re- 
spectively. The binary fraction of extremely low-mass 
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Fig. 4. — The observed binary fractions of WDs as a function 
of mass, compared to a theoretical binary evolution model (Nele- 
mans 2010, private communication). The 0.2 Mq WDs are from 
I Brown ct al. (2013), th e 0.4 Mp, WDs are fr om this paper, and the 
0.6 Mq WDs are from lHolberg et al.l 12008T ). 

WDs with < 0.25 M^, on the ot her hand, is 100% 
(jBrown et alj|20lO iKilic et al.ll2011l) . Hence, the binary 
fraction of WDs rises as the WD mass falls. In principle, 
we can use this relation between binary fraction and WD 
mass to constrain binary e volution models. 

iNelemans &: Tout! (|2005l) perform population synthe- 
sis calculations for single and double WD systems. The 
predicted WD binary fraction depends on at least three 
factors: the physics of common envelope phase evolution 
and the initial mass distribution and orbital period dis- 
tribution of the progenitor binaries. In these models, for 
example, a portion of the single lo w-mass WD pop ulation 
forms from He+He WD mergers (|Nelemansll2010[ ). 

Figure |4] compares the result of the Nelemans et al. 
population synthesis models with the observations. The 
models predict binary fractions of 100% and 75% for 0.2 
and 0.4 M© WDs, respectively. The observed binary 
fractions agree remarkably well with model predictions 
for low-mass WDs. Given the uncertainties in the pop- 
ulation synthesis calculations, however, it is worth com- 
paring the observed binary fraction of low-mass WDs 
with that of single subdwarf stars. 

The majority (90%) of Hc+He WD mergers are ex- 
pected to create subdwarf B stars, the mass of which are 
expected to exceed the helium ignition limit of ks 0.45 
Mq. The bi rth rate of single subdwarf B stars is 2x 10~ 14 
pc -3 yr _1 (Nelemans. 2010). In comparison, the forma- 
tion rate of 0.4 M w WDs in the P G survey is ~ 4 x 10~ 14 
pc -3 yr~ _1 (<Liebert et al.l [2005). Thus the formation 
rate of single low-mass WDs, given our observations, is 
i~ 10~ 14 pc -3 yr _1 . This is comparable to the birth rate 
of single subdwarf B stars, which means that the merger 
scenario by itself cannot explain the formation rate of 
si ngle low-mass WDs. 

INelemans &: Taurisl ((1998) discuss the possibility of 



creating single low-mass WDs through common enve- 
lope evolution of solar mass stars with massive plan- 
ets or brown dwarfs in close (< 3 yrs) orbits. However, 
the absence of stellar or massive brown dwarf (> 40Afj) 
companions to half a dozen apparently single low-mass 
WDs restricts the potential companions to a narrow 
mass range, where they are massive enough to expel the 
stellar envelope during the common envelope phase but 
also small enough to avoid detection in the mid-infrared 
(jKilic et all 12010) . A well tuned scenario involving a 
common envelope phase between a 1 M© star and a 20- 
40 Mj brown dwarf at a certain orbital separation can 
explain the single low-mass WDs, but it seems unlikely 
t hat this scenario w ill explain all single low-mass WDs. 

IKilic et al.1 (2007) propose that a significant fraction of 
nearby field stars have super-solar metallicity and that 
the single low-mass WD population can form through 
enhanced mass loss from these stars. We find up to 9 
single low-mass WDs fro m the 348 DA WD s in the PG 
survey, or 2.6%. Similarly. INelemans! (|2010f ) find that 15 
of 636 WDs (2.4%) observed in the S PY survey are sin- 
gle low-mass WDs. IKilic et all (|2007l ) estimate that the 
fraction of metal-rich stars with [Fe/H] > +0.3 has been 
more than 2-3% in the past 10 Gyr. Thus the observed 
frequency of single low-mass WDs is consistent with the 
enhanced mass-loss scenario. 

While binary frequency provides a useful constraint 
on formation scenarios, the period distribution is funda- 
mental to evolution of the low-mass WD binary systems. 
Figure [5] compares the observed orbital period distribu- 
tion of two complete, magnitude-limited samples: the 
0.4 Mq WDs s tudied here, and th e 0.2 M© WD sys- 
tems studied by I Brown et al.l (|2010D . Interestingly, none 
of the 0.2 M© WD systems have orbital periods longer 
than a day, whereas at least 40% the 0.4 M© WD sys- 
tems have orbital periods longer than a day. The tightest 
main-sequence binary systems will interact the earliest in 
their evolution; they will go through common envelope 
phases, lose their envelopes before helium ignition in the 
core, and end up as He-core WDs in even tighter binary 
systems. Thus finding the lowest mass WD binaries in 
the most compact binary systems is consistent with our 
understanding of binary evolution. 

4.2. Implications 

We close by noting that old, metal-rich stellar pop- 
ulations dominate in many astrophysical environments, 
namely, the bulges and spheroids of galaxies. The ex- 
istence of single low-mass WDs in the field, as well as 
the observed stellar population of NGC 6791, demon- 
strate that old, metal-rich stars are evolving through at 
least one of the unique channels discussed here. These 
evolution channels will affect the red giant luminos- 
ity function, blue horizontal branch morphology, and 
the integrated ligh t of the overall s tellar population 
(jKalirai et al.ll2007l) . lHan et al.l (|20l0) . for example, ar- 
gue that hot subdwarf stars formed from He-He WD 
mergers naturally explain the UV-upturn seen in the in- 
tegrated light of elliptical galaxies. Thus the existence of 
single low-mass WDs in the field, and the evolutionary 
channels they are linked to, has important implications 
for interpreting the age, metallicity, and mass-to-light ra- 
tio of metal-rich stellar populations throughout the Uni- 
verse. 
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Fig. 5. — The period distribution of our magnitude-limited sam- 
ple of 0.4 Mq WDs (solid histogram) compared to the ma gnitude- 
limited sample of 0.2 M WDs from lBrown et al.l <201Ch (dashed 
histogram). 



5. CONCLUSIONS 

We discuss radial velocity and photometric observa- 
tions of 30 low-mass WDs identified in the magnitude 
limited PG survey. We identify a total of 21 binary 
systems and provide orbital parameters and limits on 
companion masses. Our sample includes two WD+WD 
merger systems and two pre-cataclysmic binaries, one of 
which is a new discovery. 

Nine objects in our survey do not show significant ra- 
dial velocity variations. Six of these objects have mid- 
infrared photometry and none show mid-infrared excess 



from stellar or massive brown dwarf companions. Thus 
the fraction of single, low-mass WDs is 20%-30%. We 
discuss the potential formation channels for these single 
low-mass WDs: binary mergers of lower-mass objects, 
mass loss due to interactions with substellar companions, 
and mass loss due to super-solar metallicity winds. Bi- 
nary mergers of two lower-mass WDs would most likely 
create He-burning sdB stars more massive than 0.4 Mq. 
The birth rates of single subdwarf B and low-mass WD 
stars are comparable. Therefore, the merger scenario is 
unlikely to explain all of the single low-mass WDs. In- 
stead, a significant fraction of the single low-mass WDs 
may form as a result of interactions on the red giant 
branch with substellar companions or because of en- 
hanced mass loss from the most metal-rich stars in the 
solar neighborhood. 

Our understanding of the different formation channels 
for single low-mass WDs will benefit from infrared obser- 
vations that can put strong limits on potential substel- 
lar companions and from theoretical studies on common 
envelope evolution and the fate of short period MS + 
brown dwarf (or massive planet) systems. The observed 
mass and period distributions of our targets are useful 
for constraining binary evolution theory via population 
synthesis studies, which can probe different evolutionary 
channels for low-mass WD formation. 
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APPENDIX 
RADIAL VELOCITY DATA 

Table 2] presents our 720 radial velocity measurements for 21 low-mass WDs. We also present 11 emission line 
velocity measurements, labeled PG 1458(Mg). The table columns are: (1) Object name, (2) heliocentric Julian date, 
and (3) heliocentric radial velocity. Table [4] is available in its entirety in machine- readable form in the online journal. 
A portion of the table is shown here for guidance regarding its form and content. 
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TABLE 2 
Binary Orbital Parameters 



Object 


P (d) 


K (km/s) 


7 (km/s) 


To (days + 2450000) F-Test 


MF (Mq) 


M 2 (M Q ) 


T (Gyr) 


IN 


PG 0132+254 




C221 


17 


No detectable period 








1 


PG 0237+242 


0.7417+0.0269 


78+9 


—14+4 


4386.6917+0.0154 7.00e - 06 


0.0371+0.0134 


> 0.25 


< 190 


2 


PG 0808+595 




21+7 


28+3 


7.16e - 02 








1 


PG 0834+501 


1.2849+0.0564 


58+9 


-13+5 


4464.0581+0.0441 1.47e - 05 


0.0262+0.0132 


> 0.22 


< 920 




PG 0846+249 




(91) 


-35 


High dispersion 








4 


PG 0934+338 


1.1142+0.0055 


111+17 


14+10 


4465.7829+0.0205 5.71e - 07 


0.1604+0.0759 


> 0.50 


< 320 




PG 0943+441 




(20) 


49 


No detectable period 








1 


PG 1022+050 a 


1.157+0.001 


75+1 


39+1 




0.0500+0.0020 


> 0.30 


< 480 




PG 1036+086 


1.3283+0.0109 


83+18 


93+11 


4474.7431+0.0499 1.55e - 04 


0.0815+0.0525 


> 0.37 


< 610 




PG 1101+364 b 


0.145+0.001 


71+2 


39+1 






0.37 


2 


4 


PG 1114+224 


0.3198+0.0147 


34+7 


44+3 


4476.7679+0.0113 1.46e - 03 


0.0014+0.0009 


> 0.07 


< 64 




PG 1202+608 c 


1.493+0.001 


77+8 


0+6 




0.0720+0.0200 


> 0.34 


< 93 




PG 1210+141 d 


0.642+0.001 


131+3 


15+2 




0.1490+0.0100 


> 0.46 


< 90 




PG 1224+309° 


0.260+0.001 


112+14 


-4+12 




0.0380+0.0040 


0.28 


< 10 


2 


PG 1229-013 




22+7 


22+4 


2.48e - 02 








1 


PG 1241— 010 f 


3.347+0.001 


68+1 


16+1 




0.1110+0.0050 


> 0.42 


< 6700 




PG 1249+160 




(16) 


1 


. . . No detectable period 








1 


PG 1252+378 




30+8 


26+4 


3.00e - 02 








2 


PG 1317+453 f 


4.872+0.001 


64+1 


—24+1 




0.1320+0.0040 


> 0.44 


< 19000 




PG 1320+645 




(110) 


-16 


High dispersion 








3 


PG 1415+133 




20+8 


62+4 


7.97e - 02 








2 


PG 1458+172 


0.1653+0.0003 


143+8 


73+12 


4911.8550+0.0002 1.08e - 04 


0.1805 


0.55 


< 1.7 


2 


PG 1458 (Mg) 


0.1653+0.0003 


106+14 


13+22 


4913.7613+0.0003 6.74e - 03 


(0.0386) 


(0.21) 


(< 5.3) 






0.8603+0.0200 


45+9 


Q1 1 A 


4270.1663+0.0636 8.90e - 04 


0.0085+0.0051 


> 0.14 


< 460 




PG 1554+262 




39+10 


0+5 


1.58e - 02 








1 


PG 1614+136 




25+6 


-7+3 


1.44e - 02 








1 


PG 1654+637 




25+9 


29+4 


4.63e - 02 








1 


PG 1713+333 f 


1.127+0.001 


56+1 


56+1 




0.0200+0.0010 


> 0.19 


< 700 




PG 2226+061 




(23) 


31 


No detectable period 








1 


PG 2257+162 




27+7 


25+3 


2.30e - 02 








2 


PG 2331+290 f 


0.166+0.001 


156+3 


-11+3 




0.0660+0.0040 


> 0.34 


< 2 




References. — 
IT995I) 


(a) IMoralcs-Ruoda et all (1200511 


, (b)|Mars 


hi 119951), fcllHolbers et all 119951), (d) INelemans et all 120051). (e) lOrosz et al.1 119991). (f ) IMarsh 



Note. — (1) Probable single WD. (2) Infrared color excess, M dwarf companion. (3) Velocity variable, unconstrained orbit. (4) Double-lined spectroscopic bii 



